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Abstract 

In this paper, we apply Osgood’s criterion from the theory of ordinary differential equa¬ 
tions to detect finite-time singularities in a spatially fiat FLRW universe in the context of a 
perfect fluid, a perfect fluid with bulk viscosity, and a Chaplygin and anti-Chaplygin gas. In 
particular, we applied Osgood’s criterion to demonstrate singularity behaviour for Type 0/big 
crunch singularities as well as Type Il/sudden singularities. We show that in each case the 
choice of initial conditions is important as a certain number of initial conditions leads to finite¬ 
time, Type 0 singularities, while other precise choices of initial conditions which depend on 
the cosmological matter parameters and the cosmological constant can avoid such a finite-time 
singularity. Osgood’s criterion provides a powerful and yet simple way of deducing the exis¬ 
tence of these singularities, and also interestingly enough, provides clues of how to eliminate 
singularities from certain cosmological models. 


1 Introduction 


The Osgood criterion [T] is a classical criterion, due to W.F. Osgood in 1898 which gives conditions 
for ordinary differential equations to admit unique solutions and also to have singularities where 
solutions explode in a finite time. Following the conventions in [2], Osgood’s criterion states a 
solution x{t) of the initial value problem 


i = f 

x{to) = U 


blows up in finite time if and only if 




ds 

f{s) 


< oo. 


( 1 ) 

( 2 ) 

(3) 


This condition can be derived from Barrow’s formula [3], which says that the solution (j){t) of 
such an ODE is given by solving 

t-to= (4) 

A f{s) 
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Suppose the solution becomes singular in a finite time, tg, that is, 
we have that: 

°° ds 

+to <oo. 

Further, a solution to Eqs. ([I])-© is unique if 



Certainly, this is a criterion that is of tremendous importance for physical applications, and as 
we show in this paper, it is of particular importance in cosmology. In cosmology, assumptions of a 
spatially homogeneous spacetime, allow Einstein’s field equations to be written as an autonomous 
system of ordinary differential equations. Divergent solutions, that is, solutions that “blow up” in 
a finite time are indicative of finite-time physical and curvature singularities that occur in universe 
models. As we show below, Osgood’s criterion provides a powerful and yet simple way of deducing 
the existence of these singularities, and also interestingly enough, provides clues of how to eliminate 
singularities from certain cosmological models. Despite this, we note that Osgood’s criterion has 
received very little attention in the physics and cosmology community. For this reason, we believe 
that the work here is unique, and could be of interest in future cosmological studies. 

In this paper, we will attempt to study the question of whether an expanding spatially fiat 
Friedmann-Lemairte-Robertson-Walker (FLRW) universe with cosmological constant and ordinary 
and exotic matter will expand forever in the future or develop a finite-time singularity. It is 
understood that the ultimate fate of the universe depends on possible decay of dark energy in the 
future, and if it does not occur, the universe will expand forever [4]. 

Further, the concept of the existence of singularities in the context of general relativity has 
been studied quite extensively. With respect to the cosmological case, there is the fundamental 
singularity theorem g [5] 0 which describes irrotational geodesic singularities, and states that if 
A<0,/i-|-3p>0, and ^ + p > 0 in a fluid flow where /i is the energy density of the fluid, p is 
the pressure of the fluid, in addition to having h = 0, a; = 0, and Hq = 0 at some time sq, then 
a spacetime singularity, where either the expansion scalar goes to zero or the shear scalar diverges 
occurs at a finite proper time tq < 1/Hq before sq. 

Further elaborations are built upon this singularity theorem. It is of interest to note that in g 
five possible routes to avoid the conclusions of this singularity theorem are discussed in detail. 
They are a positive cosmological constant, acceleration, vorticity, an energy condition violation, 
or alternative gravitational equations. We refer the interested reader to g for further discussions 
regarding these issues. 

Following g, we also note that singularities occur in cosmology not only in the context of 
FLRW models, but also for realistic anisotropic and inhomogeneous models of the universe in 
which the strong energy condition /r-|-3p > 0 is satisfied. Related to this, Penrose g did pioneering 
work on black hole singularities producing a theorem that allowed one to predict the existence 
of singularities in realistic gravitational collapse cases. Hawking [8] extended these results to the 
cosmological context leading to the famous Hawking-Penrose singularity theorem which implied 
that space-time singularities are to be expected if either the universe is spatially closed or there is 
an ‘object’ undergoing relativistic gravitational collapse (existence of a trapped surface) or there 
is a point p whose past null cone encounters sufficient matter that the divergence of the null rays 
through p changes sign somewhere to the past of p (i.e. there is a minimum apparent solid angle. 



= oo, where tg < oo, then 


(5) 
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as viewed from p for small objects of given size). The theorem applies if the following four physical 
assumptions are made: (i) Einstein’s equations hold (with zero or negative cosmological constant), 
(ii) the energy density is nowhere less than minus each principal pressure nor less than minus the 
sum of the three principal pressures (the ‘energy condition’), (iii) there are no closed timelike curves, 
(iv) every timelike or null geodesic enters a region where the curvature is not specially aligned with 
the geodesic. (This last condition would hold in any sufficiently general physically realistic model.) 

Further, Bekenstein studied exact solutions of Einstein-conformal scalar equations and pre¬ 
sented a class of FLRW models which contained both incoherent radiation and a homogeneous 
conformal scalar field which bounce and never pass through a singular state, thereby circumventing 
the singularity theorems by violating the energy condition. Parker and Fulling m considered a 
classical gravitational field minimally coupled to a quantized neutral scalar field possessing mass. 
They concluded that quantum effects can sometimes lead to the avoidance of the cosmological sin¬ 
gularity, at least on the time scale of one Friedmann expansion. Collins and Ellis m examined in 
detail the singularities that occur in Bianchi cosmologies. Barrow and Matzner [12] showed that the 
singularity corresponding to homogeneous and isotropic universes observationally equivalent to ours 
must be of simultaneous Robertson-Walker type containing only small curvature fluctuations. In an 
interesting application of Robertson-Walker singularities, Barrow m showed that the large-scale 
velocity and vorticity fields required for the vortex and spinning-core theories of galaxy formation 
can be generated primordially in a natural way from a FLRW singularity. 

Related to future finite-time singularities that develop in FLRW models is the question of the 
so-called “ultimate fate of the universe”, there are several possibilities [T3] which we list in Tabled] 


Singularity Type 

Description 

T 

K 

Big Crunch (Type 0) 

9,p,p ^ 00, 1 —>■ 0 

Strong 

Strong 

Big Rip (Type I) [15] [16] 

l,p, p ^ oo 

Strong 

Strong 

Sudden Singularities (Type II) |17H37j 

1,9, p < oo,p —>■ oo 

Weak 

Weak 

Big Freeze (Type III) [38] 

1 < oo, 9, p,p ^ oo 

Weak 

Strong 

Generalized Sudden / Big Separation (Type IV) [39] 

l,9,p,p<oo,p^oo 

Weak 

Weak 

w-Singularities (Type V) |40H42| 

l<oo,p,p^0,w^oo 

Weak 

Weak 


Table 1: A classification of known finite-time singularities. I refers to the scale factor in the FLRW 
metric, p and p refer to the pressure and energy density of the matter content in the FLRW 
spacetime, w is the barotropic index, and 9 is the expansion scalar as defined above. Note that, 
sudden singularities are weak singularities where the spacetime can be extended after the singular 
event. The singularities are classified as being strong or weak according to the classification due to 
Tipler (T) [43] and Krolak (K) [44] . 

Further, so-called infinite-time singularities may also occur as well. These include Type oo or 
directional singularities, which are of strong type according to both [33] and [33]. There are also 
other types of non-singular future behaviour which are also of the infinite-time variety. These 
include the little rip [351131], pseudo-rip [33, and the little sibling of the big rip |3H]- These are 
discussed in more detail in [39] and [14] . 

It is also important to note that conditions describing the avoidance of finite-time singularities 
have been given in [18], [501153] . 

In this work, we will consider a spatially homogenous and isotropic k = 0 FLRW universe with 
a cosmological constant and different types of matter sources. We will use Osgood’s criterion to 
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establish when such models admit hnite-time, Type 0 singularities as discussed below. 


2 The Dynamical Equations 


It is well-known that covariant derivative of the four-velocity vector Ua can be decomposed as [54j 

Ua;b ■^^^06 “1“ ^ab H” ^ab ^a^bi ( 7 ) 

where Ua is the fluid four-velocity, 9 is the expansion scalar, (Tab is the shear tensor, ujab is the 
vorticity tensor, and hab is the standard projection tensor. For spatially homogeneous universes, we 
can choose the four-velocity vector to be orthogonal to the space-like surfaces in the 3 -I- 1-spacetime 
decomposition. With this, the shear tensor aab, the vorticity tensor ojab, and acceleration iia all 
vanish. One then takes projected, symmetric, and trace-free components of the decomposition 
equation above in combination with the Einstein held equations 


^ah 2^^9ah “f ^9ab ^'^ab: 


( 8 ) 


to obtain the Raychaudhuri equation 

9 + -9^ + ■^ (b- + 3p) — A = 0, (9) 

the energy density evolution equation 


fi + 9{fi+p) = 0, (10) 

and the generalized Friedmann constraint equation 

le‘^ = K^ + A, ( 11 ) 

where p denotes the matter energy density and p denotes the matter pressure. 

In fact, we have that by the Friedmann equation dm, 

lim 9t = oo ^ lim -0^ = oo => lim \Ku.{t) -|- Al = oo => lim u{t) = oo. (12) 

Therefore, the Osgood criterion provides a way of testing for a singularity where 9 and p go to 
inhnity in a hnite time. This type of singularity is a big crunch or Type 0 singularity. In the context 
of singularity classihcation, the Type 0 singularity is considered a strong singularity [3311331113] ■ 
The advantage of considering spatially hat models is that the dynamical equations (0), (Uni) can 
be reduced to a single dynamical equation via Eq. (HU. We will therefore consider the dynamics 
of such models to be represented by an evolution equation for 0, although, the analysis in what 
follows is equally as valid for an evolution equation for p. Substituting Eq. HU) into Eq. , we 
obtain 

0 = _i(l + ^)(02_3A), (13) 

where w is the equation of state parameter. 
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3 The Existence of Finite-Time Singularities 


With the main evolution equation da in hand, we now demonstrate the existence of hnite-time 
singularities. To accomplish this, we will make use of Osgood’s criterion as defined earlier. 

To apply Osgood’s criterion to our problem, let us consider the initial value problem 


0 = -l{l + w){9^-3A), 

0 (to) = Oq. 


(14) 

(15) 


Following the introductory section and evaluating the integral in ([2]), we obtain 

2ds 

(1 + w) (s^ — 3A) 

1 _ eo-VsVA 

■\/3(l + w)V~A Oq + 



(16) 


Note that ts ^ oo only if Oq = ±-\/3-\/A, where the positive value is precisely the value of the de 
Sitter expansion. Therefore, carefully choosing the initial condition to be equal to the rate of the 
de Sitter expansion avoids a finite-time singularity where 0 —^ oo. 

The value of Eq. dT51) is finite and real if either 

{A> o}n{do > VsvT}. (17) 


Hence, for choices of Oq that satisfy these inequalities, the integral in Eq. (HU) is finite, and hence 
d —>■ oo in a finite time. 

For the sake of completeness, let us consider the case where A < 0. In this case, the integral 
Eq. (fTCl) is evaluated to be 


\/3\/A(w + 1) 


(18) 


Eor A < 0, this integral is finite and real if and only if do = 0. Therefore, for this choice of 
initial condition, this universe model will exhibit finite-time, type 0 singularities. Eurther, the 
integral diverges if do = A\/3\/A, which, as mentioned above the positive root is that of the de 
Sitter expansion. Therefore, a finite-time, type 0 singularity can be avoided by choosing the initial 
condition to be that of the de Sitter expansion. 


4 FLRW Models With Bulk Viscosity 

We now briefly consider k = 0 ELRW models with bulk viscosity. Bulk viscous models are typically 
used in early-universe cosmological models. The role of the bulk viscosity is to simply add a 
pressure term in the perfect-fluid energy-momentum tensor. The perfect-fluid energy momentum 
tensor which is imposed by the symmetry of FLRW models for fluid density pressure p, and 
four-velocity vector Ua is given by 


Tab = {p+p) UaUb + gabP- 


(19) 
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By adding a bulk viscosity term, this energy-momentum tensor is now [56] 

Tab = {^J■+p)UaUb +gabP- ^0hab, (20) 

where hab is the projection tensor, and ^ is the bulk viscosity, and p = wp. 

Following we assume an equation of state 

C = (21) 

where m is a constant and is a positive constant. 

With this definition in mind, let us define the total effective pressure as 

p = p-^0, (22) 


where p is defined as p = wp. Then, by simply substituting this definition of p into Eq. m , we 
see that Eq. becomes 




6»2-3A 


— 0^{w + 1) + 3A(w -I-1) 


(23) 


It is not possible to apply Osgood’s criterion in a general sense to Eq. (l23l) . that is, for an 
arbitrary m in Eq. p3p . We will therefore consider only specific values for m in our analysis. In 
fact, the choice of m is very important and can affect the overall dynamics of general FLRW models. 
This point was extensively studied in [57H62) . We will only consider the cases m = 0 and m = 1 in 
what follows as these are physically relevant according to m and [52]. 


4.1 Case: m = 0 

This case corresponds to the case of a constant bulk viscosity. The Raychaudhuri equation 
becomes 

■ 6*^ 30k^q 3A 0'^w 3Aw 

0 = -^ + ^5-. 

^ ^ ^ ^ ^ 

where we will use 0(to) = ^0 as an initial condition. 

Then, Osgood’s criterion becomes 


ts to — 


f 

J9o 


2ds 


(1 -I- w) — 3 (1 -I- lu) A — 3sk^ 
Evaluating this integral, we find that it is equal to 

V3y/3k^i^+4A{w+l)2 


2 log 


ts to — 


This integral only diverges for 


\/3v^3F5f+4A(71+TF 

3fc5o-2eo(™+l) 


1 


- log 


2eo(u;-l-l)-3feCo 


+ 1 


\/3y/3k^^Q + 4A(ri; -|- 1)^ 


(24) 


(25) 


(26) 


03 = 


3k^o i '/iyj 4(1 -I- wYA + 3k^^| 


2{l + w) 


(27) 
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Therefore, only for the initial condition as given in Eq. (HZl), does a unique solution exist to Eq. 
(|24l) where there is no finite-time singularity where 0 —>■ oo. Finite-time singularities therefore exist, 
for initial conditions 9o (assuming — 1 < w < 1, > 0, A > 0) if 


^ ^ 1 / -h 12A -h 12Aw2 + 2AAw , 3kCo 

° ^ 2Y (n; + l)2 ^ 2(zc + l)' 


(28) 


4.2 Case: m = 1 

The Raychaudhuri equation (l23l) becomes 

0 = i (6*2-3A) (-l-w-be^o), -l<w<l, A>0. (29) 


Specifically, we now wish to apply Osgood’s criterion to the initial value problem 

1 


e = -{0^-3A)i-l-w + 0^o), 


(30) 

Oih) = 00. (31) 

Osgood’s criterion then requires us to check the convergence/divergence behaviour of 

2ds 


Jbo “ 3A) (s^o -w-1)' 

Evaluating this integral, we find that it is equal to 


(32) 


ts — to — — 


2\/3{w + l)coth ^ (tI^) - SVA^o (21og|-5o|+log |do “ 3A| - 21og|-^o6'o -b u; -b 1|) 


3\/A (—3A^q -b -b 2 w -b 1) 
The integral in (I32p has finite real values if 




:}n{( 


or 


, 1 -b 2'u; -b I /„ /» i 

A > -^ n |6»o > V3VA1 


34" 


Now, the integral in (1321) diverges if 


00 e |V34 a| U |±Vl-b3A| U 


1 -b w 

4 


(34) 

(35) 

(36) 


Therefore, as this calculation shows, the universe model avoids a singularity where 0 —>■ oo in finite 
time under the choice of any of these initial conditions, since with these choices the integral in Eq. 
(15^ diverges. 


(33) 
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For the sake of completeness, we will also evaluate the integral in (15^ for the case of a negative 
cosmological constant, A < 0. In fact, let us define A = —A^, where A 0 S K. We then find that 
the integral in Eq. (1321) has the value 


ts 


— to 

_ N/37rMr+i) ^2 log|-Co|-2 log|w - ^ 0^*0 + 1|+ log jSA^ + 6ig|) + 2Vi{w + 1) tan"i (:^) 

3A (3A2^2 ^^2^2w + l) ■ 

(37) 


This has finite real values if 


Therefore, for these values of 9q^ the universe model will admit finite-time singularities. 
Further, the integral in Eq. will diverge if Oq is chosen such that 

1 + w 

For this choice of 0o, the universe model will not admit any finite-time singularities. 


(38) 


(39) 


5 FLRW Models With a Chaplygin Gas 


In this section, we attempt to use Osgood’s criterion to determine the existence of finite-time, big 
crunch singularities in FLRW models with what is known as a Chaplygin gas. Unlike barotropic 
matter, the Chaplygin gas satisfies an exotic equation of state [4], 

p=-^, AgR, 0 < a < 1, (40) 

where A < 0 describes an anti-Chaplygin gas, and A > 0 describes a Chaplygin gas. Finite-time 
singularities as described above have been studied a number of times in the literature in the context 
of Chaplygin gas universes [HSIHTIIBHHST] . 

Using Eq. (HOj) in Eqs. m and (HH, we obtain the following version of Raychaudhuri’s equation: 



3 



6»2-3A 



For simplicity, we will consider the case a = 1, thereby, obtaining 

■ (6»2-3A)^-9Ak2 

^ ^ 2 (02 - 3A) ’ 

9{to) = 9o. 


Therefore, applying Osgood’s criterion, we are to test the divergence of 



2 (s2 - 3A) 
9A«;2 - (s2 _ 3A)2 
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ds. 


(41) 


(42) 

(43) 


(44) 











Evaluating this integral, we get 


is to — ‘2^ \ 


tan 


-.1- 




tanh 




^0 


y V3 


Z^/Ak, — 3A 1/3 V y/An — A 


\/3\/ \/A/c + A 


( 45 ) 


Real and finite solutions corresponding to Eq. (l45]) for the case A > 0, that is for the Chaplygin 
gas, occur for A > 0 and 

\/Ak > a, 0 < do < VAk + A. (46) 

Einally, divergent solutions corresponding to Eq. (I45I) occur when 


do = '/Ak + A. 


(47) 


As can be seen from the above calculations, finite-time singularities where d —>■ 00 occur for the 
conditions as described in Eq. (l46ll . while such a singularity is avoided if do is chosen such that 
Eq. (H71) is satisfied for A > 0. Therefore, the finite-time singularity is avoided for a Chaplygin gas 
with such an initial condition. 

Considering the case of an anti-Chaplygin gas, where A < 0, we will apply Osgood’s criterion 
by setting A = —a^, where a 7 ^ 0 G K. One can see looking at the expression in Eq. (I45|) . no real 
solutions exist irrespective of the choice of principal function. Therefore, we are unable to determine 
singularity behaviour for an anti-Chaplygin gas. The physical interpretation of this result is perhaps 
that the expansion never diverges in this case. 

A further interesting case occurs in the case of a Chaplygin gas, where A > 0, but with a 
vanishing cosmological constant, A = 0. In this case, applying Osgood’s criterion to Eq. (l4^ . we 
obtain: 


is io — 




•\/3 \/~A\/k 


\f%-\/~A\fK 




(48) 


As can be seen Eq. (l4^ has finite and real values if 


0 < do < \fZ\j'/Ah 


(49) 


Therefore, for these values of do, there exist finite-time singularities. Finite-time singularities for a 
Chaplygin gas with vanishing cosmological constant were demonstrated in |88] . 

Further, finite-time singularities can be avoided if do is chosen such that 


do G {0} U 



(50) 


6 On Sudden Singularities 

We now given an example of a situation where Osgood’s criterion can be used to prove the existence 
of Type II / Sudden Singularities as described in Tabled] above. We will consider an anti-Chaplygin 
gas, with a = I as per Eq. (HO]) . We will further assume that the cosmological constant vanishes, 
A = 0. In the case of a sudden singularity, the pressure diverges, that is, p —> 00 , while the energy 
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density and expansion vanish. To apply Osgood’s criterion to this situation, we use Eqs. (UnD, (ED, 
and (1401) to formulate a new initial value problem for the pressure, p{t): 


P = 


2 \ - 1/2 
an ' 


VSk (a^ + 


p{to) = po, 


(51) 

(52) 


where for simplicity, we have denoted the anti-Chaplygin gas constant A, hy A = —a?', where 
a ^ 0 G ffi.. 

Applying Osgood’s criterion to the initial value problem (l5T]) - (l54D . we are to evaluate 


ts — to — 




1/2 


-ds. 


(53) 


loo \/3 (a'^K + KS^) 

We find that 

ts-to= (54) 

^/d - log\y/2^/dy/p^ -I- a -|-po| + log \-\/2^y/p^ -b a -|-po| - 2tan“^ + l) + 2tan“^ ^1 


2\/6ay/K 


■log|3-2V2| 


One can easily confirm that Eq. (l54l) is real and finite for all po G M. Therefore, for all po S K, 
we have that the universe model will exhibit a finite-time sudden singularity. 


7 Conclusions 

In this paper, we have applied Osgood’s criterion to detect finite-time singularities in a spatially 
flat FLRW universe in the context of a perfect fluid, a perfect fluid with bulk viscosity, and a 
Chaplygin and anti-Chaplygin gas. In particular, we applied Osgood’s criterion to demonstrate 
singularity behaviour for Type 0/big crunch singularities as well as Type Il/sudden singularities. 
We have shown that in each case the choice of initial conditions is important as a certain number 
of initial conditions leads to finite-time singularities, while other precise choices of initial conditions 
which depend on the cosmological matter parameters and the cosmological constant can avoid a 
finite-time singularity. 

As we have shown, Osgood’s criterion provides a powerful and yet simple way of deducing the 
existence of these singularities, and also interestingly enough, provides clues of how to eliminate 
singularities from certain cosmological models. For this reason, we believe that the work here is 
unique, and could be of interest in future cosmological studies. 


V2y/po\ 
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